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Abstract 

We investigate the effective case of two-fiavor neutrino oscillations in infinitely dense 
matter by using a perturbative approach. We begin by briefly summarizing the con- 
ditions for the three-flavor neutrino oscillation probabilities to take on the same form 
as the corresponding two-flavor probabilities. Then, we proceed with the infinitely 
dense matter calculations. Finally, we study the validity of the approximation of 
infinitely dense matter when the effective matter potential is large, but not infinite, 
this is done by using both analytic and numeric methods. 
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1 Introduction 



The large experimental collaborations [1-10] in neutrino oscillation physics 
have successfully used effective two-flavor neutrino oscillation formulas in their 
analyzes in order to determine the fundamental neutrino oscillation parame- 
ters. However, we now know that there exists (at least) three active neutrino 
flavors in Nature, which, in principle, means that the analyzes have to be 
carried out using three-flavor neutrino oscillation formulas. Thus, it is impor- 
tant to carefully examine the effective two-flavor formulas in order to deter- 
mine their validity in different situations, especially since neutrino oscillation 
physics has now entered the era of precision measurements. For example, the 
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Super-Kamiokande [1-3], SNO [4-6], K2K [7], and KamLAND [8-10] collabo- 
rations have significantly pinned down the errors on the fundamental param- 
eters and future results of the above mentioned experiments as well as other 
long-basehne experiments will continue to decrease the uncertainties of these 
parameters. 

Furthermore, it is known that the presence of matter can result in large al- 
terations in the behavior of neutrino oscillations [11]. An example of this is 
the Mikheyev-Smirnov-Wolfenstein (MSW) effect [11,12], which is the most 
plausible description of the oscillations of solar electron neutrinos into neu- 
trinos of different flavors. The matter effects are proportional to the energy 
of the neutrinos as well as the matter density. Thus, for large neutrino en- 
ergy and dense matter, the matter effects become important. Note that, in 
Ref. [13], approximate mappings between the three-flavor neutrino parameters 
in vacuum and matter have been obtained. In addition, in Refs. [14-21], ex- 
act treatments of three-flavor neutrino oscillations in constant matter density 
have been performed, whereas, in Refs. [13,22-34], different analytic approx- 
imations of three-ffavor neutrino oscillation probability formulas have been 
investigated for both constant and varying matter density. 

In this Letter, we study the matter effects on neutrino mixing and oscilla- 
tions in the limit where the matter density becomes infinite and the resulting 
neutrino oscillation probabilities are actual two-fiavor formulas. We also in- 
vestigate how well these two-fiavor formulas reproduce the actual neutrino 
oscillation probabilities when neutrinos propagate through matter with large, 
but not infinite, matter density. 



2 Effective Two-Flavor Formulas 



With three neutrino flavors, there are six fundamental parameters which in- 
fluence neutrino oscillations, two mass squared differences {Arri^i and Amli), 
three mixing angles {612, 623, and ^^13), and one CP- violating phase (6). For 
the leptonic mixing matrix U, we adopt the standard parameterization given 
in Ref. [35]. In addition, we introduce the mass hierarchy parameter a = 
Am|i/Am|i as well as the parameter A = Am|i/(2£'), where E is the 
neutrino energy. We also use the abbreviations Sij = sin^j^, Cij = cos 9ij, 
ss = sin S, and cs = cos 5. 

In some special cases, it is possible to write one or more of the neutrino 
oscillation probabilities as effective two-flavor formulas, i.e., the probabilities 
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can be written as 



Pap = <^a/3 + (1 - sm'(2^) sin^ ( -^L\ , (1) 



where P^p is the probabihty of the transition ~^ ^0, sin^(26') defines the os- 
cillation amphtude, Am^ defines the oscillation frequency, and L is the length 
of the path travelled by the neutrinos. In the case of three-fiavor neutrino 
oscillations in vacuum, the survival probability Paa is given by 

Paa = 1-4 E \U^^\'\UaJ\'sin'(^L]. (2) 



l<i<j<3 



This probability takes the form of Eq. (1) if any of the elements Uai of the 
leptonic mixing matrix or any of the mass squared differences equals zero. 
The same argument holds for the neutrino oscillation probability P^^. If one 
of the conditions is satisfied, then the leptonic mixing matrix can be made 
real and there is no CP-violation in neutrino oscillations (Majorana phases, 
which do not influence neutrino oscillations, could still be present if neutrinos 
are Majorana particles). In that case, the neutrino oscillation probability is 
given by 



Amf, 



Note that if Uai — 0, then only the oscillations involving the neutrino flavor 
will become effective two-flavor cases. 



The two-flavor cases mentioned above are realized if we make either of the 
approximations 9 is or a — > 0, see Table 1 for amplitudes and fre- 
quencies. These approximations can be used when analyzing neutrino os- 
cillation experiments neglecting matter effects, for example, the oscillations 
of electron anti-neutrinos of reactor experiments such as KamLAND [8-10] 
(^13 0) and oscillations of atmospheric neutrinos in experiments such as 
Super-Kamiokande [1-3] (a — > 0). 

The above argument also holds for neutrinos propagating through matter of 
constant density if exchanging the fundamental neutrino parameters for their 

effective counterparts in matter. This gives us another important two-flavor 
case, namely the limit of an effective matter potential approaching inflnity, 
where all neutrino oscillation probabilities reduce to the two-flavor form of 
Eq. (1). The remainder of this Letter is devoted to study this hmit in detail. 
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Tabic 1 

The neutrino oscillation amplitudes and frequencies for different two-flavor formulas 
in vacuum. In the a = case, all channels are two-flavor channels and the amplitudes 
and frequencies can be deduced from the channels in this table (see Ref. [34]). 

3 Mixing in Dense Matter 



In matter, the Hamiltonian of three flavor neutrino evolution consists of two 
parts, the kinematic term H]^ = AUdiag{0,a, 1)W and the interaction term 
Hjn — diag(y, 0,0), where V — \/2Gpne, Gp is the Fermi coupling con- 
stant, and Ue is the electron number density, resulting from coherent forward- 
scattering of neutrinos. When 2EV ^ Amg^, the term Hk can be thought of as 
a perturbation to the term Hm- Since Hm has two degenerate eigenvectors, we 
can use degenerate perturbation theory to find approximate eigenvectors and 
eigenvalues of the total Hamiltonian. The degenerate sector of Hj^ is spanned 
by — (0, 1, 0)^ and u-,- — (0, 0, 1)^. In this sector, the Hamiltonian is of the 
form 

, a b* , 

H,= \ 1, (4) 
b d 

where a, b, and d depend on the fundamental neutrino parameters as well as 
the effective matter potential V and the neutrino energy E. 

In the case of two-flavor neutrino oscillations in vacuum, the Hamiltonian 
takes the form 



H 



Am2 / - cos(2e) sin(2^) ^ 



4^ i sin(2^) cos(2^) 



(5) 



where Am^ and 6 are the mass squared difference and mixing angle in the 
two-flavor scenario, respectively. Defining 6^ as the effective mixing angle in 
the degenerate sector of if^ and Am^ as the effective mass squared difference. 
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Fig. 1. The quotient |Am^/Arn32| as a function of the leptonic mixing angle ^13 for 
different values of the mass hierarchy parameter a. In this plot, the mixing angle 
9i2 has been set to 33°. 

by comparing Eqs. (4) and (5), we obtain 



sin2(2^^) 



4|6P 



d)2 + 4|6|^ 



= 4 



{523C23[Ci3 - ^(€^2 - -§12^13)] ~ «-Sl2Cl2Si3Q COs(26'23)}^ 
[Cl3 - a(Ci2 - Sl3Sl2)]^ + 4ac?2S?2Si3 



+4 



'^^'^5'^12Cl2'5l3 



13 



«(Cl2 - S?3sf2) 



+ 4:acl2Si2Si3 ' 



Am^ = 2EJ(a-rf)2 + 4|6|2 



= I Am^i I ^[cfg - a(ci2 - sfasfs)]^ + 4Q;c?2sf2S?3. 



(6) 



(7) 



We note that even for 6'i3 = or a = 0, the effective two flavor parameters 
Am^ and 9oo are not equal to Am32 and 6*23, which might have been expected. 
However, for reasonable values of the fundamental neutrino parameters, ^00 
does not deviate significantly from ^23- In Fig- 1, the quotient lAm^/Amlgl 
is plotted as a function of 613 for different a. From this figure, we conclude 
that the sign of the correction to the approximation Am^ = | Am^gl depends 
on the values of the parameters 613 and a. It should be noted that if the ab- 
solute value of the fundamental parameter Am|2 and the effective parameter 



AmL were to be measured 



I.e.. 



Am^ determined in an experiment where 



the matter is extremely dense and |Am32| is determined by vacuum or close 
to vacuum experiments), then the quotient |Am^/Am32| could provide valu- 
able information on the mass hierarchy {i.e., the sign of a) and the leptonic 
mixing angle ^13. It is also interesting to observe that for a > (normal mass 
hierarchy) , there exists a value of ^13 such that Am^ 



|Ami2|. 



In Fig. 2, we plot the effective neutrino parameters in matter as a function 
of the product EV. For neutrinos at large values of EV, the effective mixing 
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Fig. 2. Numerical results for the effective neutrino parameters in matter as a function 
of EV for a normal mass hierarchy with the vacuum parameters Am^i = 8 • 10~^ 
eV^, Am|i = 2 • IQ-^ eV^, 9i2 = 33°, 023 = 45°, ^13 = 10°, and (5 = 0. 

angles $12 and 623, where the tilde denotes that these are the effective matter 
parameters, do not play a very important role, since 613 — > 90° when EV 00. 
Instead, in this case, the important mixing angle is the mixing between the 
first and second eigenstates, which for ^^13 = 90° becomes dependent on the 
the effective mixing angles 6^12 and 623, as well as the effective phase 5, this 
angle is clearly the angle ^00 ■ 

As can be seen from Fig. 2, we may consider essentially three different regions 
for the parameter EV. The first region is 2EV <C Am^i, where the funda- 
mental neutrino parameters are essentially equal to the vacuum parameters, 
the second region is 2EV ^ Amg^^, where is essentially an eigenstate to the 
Hamiltonian and there is two-flavor neutrino mixing between z/^ and Ur- The 
last region is the region inbetween the other two, where resonance phenomena 
occur. In general, the results for an inverted mass hierarchy are similar to 
those of a normal mass hierarchy. 

Observe that for atmospheric neutrinos passing through the Earth with a char- 
acteristic energy of ii^ ~ 1 GeV, the product EV is typically in the upper part 
of the resonance region in which none of the effective mass squared differences, 
nor any of the effective mixing angles, are small. This effect, which is studied 
in, for example, Ref. [36] (for a recent study, see Ref. [37]), should be taken 
into account when analyzing data from atmospheric neutrino experiments. 



4 Accuracy of the Two-Flavor Approximation 



The accuracy of the two-flavor approximation in the limit EV — > cxo is depen- 
dent on how close the leptonic matter mixing angle ^13 is to 90°, or in other 
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words, the value of C13. The correction to U3 — Ue is given by 



1^3 -l^e- 



V 



(8) 



in first order non-degenerate perturbation theory (which can be used, since i/e 
is a non-degenerate eigenstate of Hm), where H^a is the e-a element of Hk. 
Prom this calculation follows that to first order in Am|i/(£'y), the quantity 
Ci3 is given by 



In Figs. 3 and 4, wc plot the relative and absolute accuracies of the large matter 
density approximation for the probability P of a z/^ to oscillate into some other 
flavor in the limit EV 00 for two different values of L/i? as a function of 
EV and 613. Since this oscillation probability is given by P = 1 — P^^, the 
relative accuracy is given by ^rei = ^abs/-P = ^abs/(l - -P^m) where ^abs is 
the absolute accuracy. As for all relative accuracies, this definition clearly has 
a problem for P = 0. Since P^^ is close to one in the regime studied in the 
plots, the relative accuracy for P^,^ is approximately the same as the absolute 
accuracy. The values used for L/E correspond roughly to values that might 
be used for future long-baseline experiments and neutrino factories. Note that 
these values are smaller than the oscillation length determined by Amg^. 

The reason why the approximation becomes worse at large values of 6*13 is 
the fact that if ^13 is large, then 6*13 will approach 90° slowly when EV — > 00 
compared to when 613 is small, see Eq. (9). Thus, larger values of EV will be 
required before the approximation EV —>■ 00 becomes valid. For EV —>■ 00, 
the region in which the two-flavor neutrino approximation Am^ = ^^^d 
6 = 623 is accurate is the region in which ^13 obtains the value for which 
Amoo — |Am32|. This is to be expected, since 9^0 — O23. When studying 
Figs. 3 and 4, one should keep in mind that the survival probability P^^^ is 
larger for L/E = 7000 km / 50 GeV than for L/E^ 3000 km / 50 GeV and 
that we have plotted the relative accuracy of the approximation. 

The reason why the two-flavor neutrino approximation EV 00 seems to be 
a good approximation even for small values of EV when ^13 is small has to do 
with the fact that for small values of L/E and 613 = 0, this approximation gives 
an error which is of the order to the {L/ E)"^ dependence of the probability 
P^fj^ when making a series expansion in the quantity L/E. 

Presently, the number of experiments in which the approximation of EV 00 
can be justified is somewhat limited. One such example could be a neutrino 
factory setup with a "magic" baseline (L ~ 7250 km) [38] with a neutrino 
energy of about 50 GeV. For this energy, the value of log(PV/eV^) typically 



Cl3 ^ 




I sin2(2^i3) - 2aslsclssl2 + a^cj^sUch + cfasL)- (9) 
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log(VE/ 1 eV^) 




log(VE / 1 ev') 



Fig. 3. The relative (upper panel) and absolute (lower panel) accuracy of us- 
ing the approximation EV ^ oo in the neutrino oscillation channel 1 — P^^ for 
L/E = 7000 km / 50 GeV. The fundamental neutrino parameters have been 
set to Am^i = 8 • lO'^ eV^, Am|i = 2 • lO'^ eV^, 6*12 = 33° 6*23 = 45°, 
and 6 = 0. The colored regions correspond to different values of the rela- 
tive accuracy ^rei = \P^l^l,nnm - P^l^l,a.pp\/{'i- - Pfiii,num) and absolute accuracy 
^abs = \PiJ,fi,num " -P^i/i,app| as indicated by the legends above the panels. The curves 
corresponds to the isocontours of the accuracy of the two flavor approximation 
9 = 023 and Am^ = Arrigj. The solid curves are the Aj-ei = 0.25 % isocontours, 
the dashed curves are the Arei = 0.5 % isocontours, and the dotted curves are 
the Aj-ei = 1 % isocontours in the relative accuracy plot. In the absolute accuracy 
plot, the solid curves are the ^abs = 0.0001 isocontours, the dashed curves are the 
^abs = 0.0005 isocontours, and the dotted curves are the 74abs = 0.001 isocontours. 

varies between —2.5 and —1.8 along the baseline, which is basically the middle 
region of Figs. 3 and 4. Even though the real Earth matter density profile is 
not constant, the EV 00 approximation should be at least as good as for 
the case of a constant matter density corresponding to the lowest density of 
the Earth matter density profile. The approximation has no direct application 
to solar or supernova neutrinos for which there is a larger V than that given 
by the Earth. This is due to the fact that these situations are well-described 
by adiabatic evolution and that the neutrino energies are lower than in the 
case of a neutrino factory. The approximation might also be applicable to de- 
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-3 -2 -1 

log(VE/ 1 eV^) 




log(VE / 1 ev') 



Fig. 4. The relative (upper panel) and absolute (lower panel) accuracy of using the 
approximation EV ^ oo in the neutrino oscillation channel 1 — P^^ for L/E = 3000 
km / 50 GeV. The areas and curves are the same as described in the figure caption 
of Fig. 3. 

termine the second order oscillation effects in high-energy neutrino absorption 
tomography of the Earth. 



5 Summary and Conclusions 

We have studied effective two-flavor neutrino mixing and oscillations in the 
limit EV —>■ oo. In this limit, the neutrino oscillations take the form of 
two-flavor ^ oscillations with the frequency and amplitude given by 
Eqs. (6) and (7), respectively. It has been noted that, if measured, the quo- 
tient lAm^/Amlal could provide valuable information on the mass hierarchy 
in the neutrino sector and the leptonic mixing angle 613. We have also ex- 
amined the validity of using the mentioned two-flavor formulas for large, but 
not infinite, matter density. The deviation from the two-flavor scenario is de- 
termined by the deviation of C13 from zero, which is given by Eq. (9) to flrst 
order in perturbation theory. The validity of the approximation has also been 
studied numerically with the results presented in Figs. 3 and 4. 
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